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Определено понятие и особенности 4ЮЇ) гравитационного мониторинга. Обработку его данньх на 
коротких профилях предложено производить по соотношениям с бьтстро убьтвающими ядрами, а ва- 
риации гравиполя из-за влияния малоинтенсивньх геофизических факторов учитьвать по методике 
Двулита. 
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1. К обоснованию мониторинга. На фоне очевидного глобального технологического про- 
гресса в условиях зкономического спада на количество, методику и возможности осуществ- 
ления крупномасштабньх геофизических работ влияют две тенденции. С одной стороньт - 
возросшая точность и производительность гравиметрических сьемок (зтот метод при прием- 
лемой точности остается недорогим поисково-разведочньтм решением в общем комплексе ге- 
офизических исследований вследствиєе совершенствования аппаратурьт и топогеодезического 
обеспечения - систем СР5). С другой стороньт заметно резкое снижение обьема измерений. 

Первое повлекло за собой необходимость пересмотра методов обработки полученньтх дан- 
ньїх, в частности, более точньг6й учет поправок Буге |Бьічков, 2007). Последнее, на фоне необ- 
ходимости обнаружения более глубинньх источников аномалий", влечет за собой пересмотр 
методики измерений с целью учета тонких особенностей гравитационньх аномалий ,щмалой 
кровью?, т.е. без усложнения математического аппарата, методик измерений и увеличения 
материально-технических затрат. Зти особенности можно ,,подсечь", если задействовать до- 
полнительную переменную - время. 

В зтом плане мировьеіе ,трендьт" геофизических наблюдений постепенно тяготеют к не- 
прерьвному 4Ї мониторингу |ВабйаєПа еї а!., 2003) исследуемой площади, изучая зволюцию 
гравитационного поля за время зсплуатации участка либо на протяжений интервала его рез- 
кой динамической активизации. Следуєт подчеркнуть, что в западной литературе | Сеорпу- 
зіс8, 2005) под термином ,,вариации сильт тяжести" имеют в виду временную разность реаль- 
ньх аномалий на ограниченньжх площадях, источниками которьгх вьтступают обьекть!т с бьіст- 
роизменяемой глубинной динамикой, в то время как в отечественной литературе за зтим по- 
нятием зарезервированьт слабьте квазипериодические флуктуации на сверхдлинньх профилях, 
пересекающих области контрастньїжх современньх вертикальньх движений блоков земной 
корь. 

Сама идея мониторинга уже витала в воздухе | Лоссовский, 2006| и некоторьг6"е ее прило- 
жения в определенной степени иллюстрирует работа |Юргин, 2006). Тем не менее, в указан- 
ньх источниках не упоминается о необходимой, по нашему мнению, подоплеке мониторинга 
- вариациях сильт тяжести в том пониманий, которое в них заложено работами Г.Т. Собака- 
ря, и их зависимости от ряда малойнтенсивньх природньх и техногенньх факторов. 

Применение повторньх наблюдений в некоторой области в гравиметрии для нас сравни- 
тельно неново (его злементь применялись в созданий региональньгїх плотностньх моделей 
Кавказа, о чем упоминаєтся в | Алексидзе, 1985, а также см. ПЮркина, 19781), хотя непрерьбів- 
ньми во времени их назвать трудно, в то время как в сейсмометрии зто давно - общеприня- 
тая практика. 

Отметим, что при организации любьтх вьшокоточньх измерений, в которьїх сила тяжести 


; Возможности методов регуляризации в решений задач построения сложньх разрезов при современном уров- 
не модельньюішх представлений о геологической среде близки к технологическому пределу. 

7 С.82 - »- ЄСЛИ СОВременньюе физико-математическиє теориий строят и исследуют свой функциональнье связи 
и последующие смьютловьюе рядью на территорий преимущественно пространства, то зксперимент зондируєт 
пока еще не подвластное теорий время". 


- отдельньй параметр, ее изменение следует обязательно исследовать и учитьвать |Соба- 
карь, 1972. 

2. О вариациях силь тяжести. В |Собакар, 1972) впервьге в отечественной литературе 
проанализированьт различнье составляющие вариативной части гравитационного поля и вь- 
делень неприливньте квазипериодические вариации (КТІВ) сильт тяжести (с аплитудой в пре- 
делах 3-5-кратной погрешности измерений). Для их измерений бьмла создана новая методоло- 
гическая (неоднородная модель Земли как в физико-химическом, так и в знергетическом пла- 
не), методическая (КПВ максимальнь при пересечений разновозрастньх тектонических 
структур и в областях контрастньх современньтх вертикальньх движений) и метрологичес- 
кая (наблюдения на сети оптимальной плотности и конфигурации несколькими приборами с 
вьшсоким козффициентом надежности и сравнение градиєнтов КПВ и наблюденньх анома- 
лий) база. 

Сами КПВ очевидно тесно связань с зндогенньми процессами формирования и разви- 
тия плотностньїх неоднородностей земной корь и мантии, хотя не исчерпьваются ими. В хо- 
де перераспределения вещества внутри Земли меняется ее момент инерции, а как следствие - 
ротационньй режим, и с ним сила тяжести. При зтом перераспределениий нарушаєтся фигура 
равновесия Земли, что изменяєт силу тяжести. В ходе изостатического восстановления фигу- 
рьт равновесия сила тяжести снова меняєтся. Зта гамма процессов и генерируєт КПВ силь 
тяжести Земли. Обратимая часть процессов создаєт периодическую часть вариаций, а необ- 
ратимая часть - непериодическую часть вариаций, которая и формируєет стационарное грави- 
тационное поле Земли (гравианомалий). 

Отмечена |Лубков, 2010) корреляция крупномасштабньтх плотностньх неоднородностей 
мантии и приливньтх параметров Земли, характеризирующих ее приливное деформированиє. 
Сами приливньюе параметрь (т.наз. числа Лява и Шида) находят на основе расчета релакса- 
ционньх амплитуд гравипотенциала на поверхности Земли и ее ядра. В связи с зтим возмож- 
но наличие кросскореляции КПВ и упомянутьх динамических параметров Земли, что требу- 
ет отдельного изучения. 

Невзирая на слабую интенсивность, аномалий КПВ сильт тяжести вьшявляют вполне уве- 
ренно вследствие особенностей поведения временной кривой КІВ и тесной обратнопропор- 
циональной корреляции с кривьгми вертикальньх движений земной корьг. Последнюю счита- 
ют |Собакарь, 1972) следствием общности процессов в верхней мантии, влияющих на КПВ и 
вертикальньєе движения. Суммарная величина КТПВ является суммой суперпозиций вариаций 
различного происхождения, знака, периода и апмлитудь. Зта суммарная амплитуда в конеч- 
ном итоге и определяет зволюцию гравитационного поля Земли, вьзванную зволюцией 
неоднородностей ее корьг и верхней мантии. 

Ввиду сказанного, вариации силь тяжести следует учитьвать в интерпретации резуль- 
татов 4Ї) гравиметриий - с целью введения соответствующих поправок в гравитационнь"е 
сьемки разньх лет, в долговременную прецизионную топосьемку, уточнения реологиий ис- 
следуемого участка и др. 

3. Основа мониторинга. Определим основнье понятия. Назовем гравитационньм мони- 
торингом ряд периодически повторяемьх в реальном времени непрервівньх на протяжений 
фиксированного временного отрезка (рис. 1) микрогравиметрических измерений гравитаци- 
онного поля и обработку его данньх с учетом влияния окружающей средьт и сферьт примене- 
ния. 

Величина временного отрезка зависит от качества измерений, мерьї неопределенности ре- 
зультатов наблюдений, динамики (амплитудьт и частотьт) поля. 

Физической основой гравиметрического мониторинга является непрерьгвная связь дина- 
мики гравитационного поля и параметров средь - ундуляциям сантиметрового диапазона от- 
вечают вариации силь тяжести в несколько мГал. Если деформация дневного рельефа опре- 
деленного участка есть прямьтм следствием приповерхностного распределения масе, грави- 
тационньй мониторинг можно применять для исследования разуплотнения и флюидного ре- 
жима зтого участка. 


Пространственное распределение вариаций вертикальной производной ИЙ, потенциала си- 
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ль тяжести прямо коррелирует с плошадньтм распределением плотностей, а временньге вари- 
ации У, четко определяют вертикальнь6е вариации насьщения флюидов. 


Аппаратная основа мониторинга средьт представляєт собой совместньєе крупномасштаб- 
ньве измерения превьшений рельефа по данньтм СР5 и абсолютньх значений силь тяжести 
(сотни стационарньх пунктов на сотни км.). Значительно дешевле обходится гравиметрия с 
относительньм учетом сильг тяжести, но ей присущи серьезнье ограничения - привязка к 
опорной сети и необходимость синхронного учета ,сползания нуля". Тем не менее, в наших 
условиях именно такие измерения при надлежащем методическом обеспечениий имеют най- 
большую перспективу. На некоторьгх участках возможно применять непрерьівнье измерения 
силь тяжести в скважинах, преобладающие в разрешающей способности наземную сьемку 
благодаря большей близости к источникам и устранению приповерхностньх влияний. Для 
компенсации влияния температурь и других внешних факторов используют некоторье на- 
работки из аппаратурной базьт морской гравиметрий. 

4. О вьщделений слабого сигнала. В целом, в исследованиях, посвященньтх 4Ю) гравимо- 
ниторингу |Сеорпубісз, 2008) , отмечена амплитуда сигнала в диапазоне от 20 мГл ГУ/1Шат5- 
опез, 2008| до 80 мГл |Сеїйпез, 2008). В границах действующих вулканов амплитуда 
сигнала возрастаєт до 300-600 мГл, а в пределах водньх бассейнов - 200-250 мГл |Фамі5 еї 
аі. 2005) при нелинейном учете ,сползания нуля". Вследствиє повторяємости измерений, 
фиксирующей невязку значений силь тяжести в точке наблюдений, можно оценить точность 
определения полезного сигнала на временном отрезке. 

С целью вьщделения слабого сигнала на фоне помех в гравиметрий применяют методьт 
корреляционного анализа и вьтчисление некоторьїхх составляющих градиентов по прямьгм из- 
мерениям силь! тяжести". Авторьг ГВайазііа еб а!., 2008) вьіделяют полезньй сигнал путем 
расчета разностной аномалий Ля, - 9, - М, - ДЕ у - ДЄ,» Где А - разность между двумя 


соседними отсчетами во времени, У « -308.6 мГал/м - поправка в свободном воздухе, /. - 
вертикальное смещение (относительное поднятие-опускание), Ля, - аномалия Буге дефор- 


мации (вклад от изменения обьема из-за сжатия средь вокруг источника, что подразумеваєт 
смещениєе плотностньх границ в неоднородной среде); Де, «212С,9б, «4296. - влияниє 


грунтовьїх вод (действие гидротермального фактора). Кроме того, используєтся привязка к 
независимьм наблюдениям на контрольном пункте, расположенном вблизи исследуемого 
участка. 

Опеределеннье улучшения слабого сигнала на фоне помех получают и путем статисти- 
ческой фильтрации влияния температурньх флуктаций с помощью алгоритмов с нечеткой 
логикой |Апдб, Сагбопе, 2004). Если не хватаєт длиньг ряда наблюдений, используют поли- 
номиальную аппроксимацию низкого порядка |Вопуаїої еї а!., 1998| с соответствующей ,,ка- 
либровкой" порядка полинома. Кроме того, используют сравнение результата фильтрации с 
данньми, полученньми с близко расположенного контрольного пункта от участка наблюде- 
ний. 

Примечательно, что как и в случае измерений КПВ силь тяжести, многие авторьт из |Се- 
орПубісє, 2008) отмечают корреляцию слабого сигнала с вертикальньтми подвижками и, пред- 
полагая разночастотность напряженности поля сильт тяжести и вариаций помех, применяют 
для усиления полезного сигнала частотную фильтрацию. Аналогию можно усмотреть опять 
же в морской гравиметрии. На равнинньх участках сигнал имеет мальй градиент, и такой 
подход, по нашему мнению, малозффективен. В таком случає нужно учитьвать влияние по- 
мех иньми методами - изменением геометриий сети наблюдений, временной фильтрацией, 
вьтчислением производньх и т.п., предлагаємьми в | Собакадь, 1972. 

5. О методике мониторинга. Вьїгскажем некоторье соображения об организации мони- 
торинга. Его методическую сторону, видимо, целесообразно основьтвать на достижениях 
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Точно так же как в морской гравимерической сьемке применяют вьічисление вертикальнькх и горизонтальньтх градиентов 
силь тяжести и их комплексную интерпретацию.Зта методика интерпретации опробована на морских нефтегазовьтх струк- 
турах |Юргин, 20061. 


| Собакарь, 19721, а изучение метрологических нюансов - на примере геодинамическкого по- 
лигона с готовой инфраструктурой и базой геофизических измерений (например, одного из 
использованньхх в исследованиях |Двуліт, 19991). Так, для мониторинга условного месторож- 
дения необходима сеть детальной рядовой сьемки масштаба 1:10000 на ровном полигоне 
10х10 км' по стандартной методике с минимумом опорньх пунктов |Гравиразведка..., 1990, 
с. 106): 21 профиль (по 20 точек с шагом 50 м) с контролем 30906 от общего количества рядо- 
вьїх точек (свьше 13 ть. пунктов измерений). Организация полигона и наблюдения на ука- 
занной сети за месяц обойдутся в -- 100 тьс. грн, без учета стоимости и аммортизации аппа- 
ратурьт и других накладньтх расходов. 

Недостаток инфраструктурь и оборудования увеличит затрать на организацию монито- 
рингового полигона на неопределенную стоимость. Использованиє цифровой регистрации 
позволит за счет упрощения схемьт сьемок обойти некоторье трудности в начальньтх зтапах 
сьемки и ускорить создание цифровой модели обьекта. Вьход из финансово-технологическо- 
го тупика - сотрудничество разньх институтов с совместньтм использованием инструмента- 
рия, персонала, методик и общим планом исследований на общем обьекте (зтим путем инте- 
грированного мониторинга (соПабогайуе дтаміту пебмогК) идут многие западньєе компаний). 

Следует отметить, что применениє классической схемь измерений на регулярной сети 
пунктов и последующий перерасчет значений по известному интегралу Пуассона пригодно 
для региональньгх исследований |РДамі58 ег а!., 2008|, но в локальньх условиях, где чаще всего 
применяют гравимониторинг |Болотнова, 2007), имеет ряд недостатков |Дубовенко, 2002). К 
тому же, иногда по ряду причин организация регулярной сети невозможна, а пересчет из не- 
регулярной сети на регулярную является задачей более сложного порядка, чем обратное вос- 
становление строения средьт по данньгтм измерений. 

Решение обратньх задач гравиметрии с данньми в псевдорегулярньгх сетях с применени- 
ем моделей средьтг типа ,бесконечньй профиль" приводит к появлению плохо обусловлен- 
ньїх систем линейньх уравнений, генерирующих бессодержательнье результать. Ввиду зто- 
го, и исходя из преимущественно краткой длиньт реальньїх профилей, целесообразен альтер- 
нативньй подход. 

Для интерпретации измерений на коротких профилях предложена |Дубовенко, 2002 сис- 
тема линейньх интегральньх уравнений с бьтстро убьтвающими ядрами: 


зма о | а ар ногі 


зб -1 
здано ЇВ нав Я анг 
п Зо п 
с) 5о (чн 5 (х)я мо), с (а) о (х)чю 55 (2), песо 
й 5 я 

С учетом зтого методика, изложенная в |Болотнова, 2007) остаєтся действенной в опреде- 
ленньх условиях (региональньй фон - полином 1-й степени; известнь плотности и положе- 
ния границ тяготеющих тел на поверхности; зти тела близки или имеют общие контакть ). 
Она, в частности, предполагает построение пространственной плотностной модели средьт в 3 
зтапа: разделение аномалий от тяготеющих тел и определение зффективньх глубин их зале- 
гания и квазиплотности (нулевое приближение), определение истинньх глубин и плотностей 
тел по решению плоской обратной задачи (первое приближение) и окончательное решение в 
среде АРО-ЗР. 

Нами предлагаєтся в методике |Болотнова, 2007) использовать программньй комплекс 
| Старостенко и др., 2004), разработанньвй на основе | Старостенко, Легостаєва, 19981, и под- 
борку программ, полученную в исследованими |Дубовенко, 2002). 

6. К интерпретации данньтх. Целью мониторинга является оценка глубинь к источнику 
аномалий и изменения обьема по данньм деформации рельефа требует знания приповерхно- 
стного распределения масе (по данньм гравиметрии). Деформации земной поверхности по- 
лучают по данньм спутниковой геодезий СР5, имеющих ряд преимуществ перед традицион- 
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ньтми геодезическими методами". Приповерхностньве неоднородности строения средь (карс- 
ть, плавунь, зонь обводнения и разуплотнения), сложная структура площади (складчатость, 
солянокупольная тектоника, разломьт), факторьт поглощения полезного сигнала (температу- 
ра, инструментальнье влияния) ограничивают зффетивность мониторинга, не снижая его 
практического значения. 

В |Ванаєна еї а!., 2008; ГДаміз еї а!., 2003) не учтена особенность вариаций гравитацион- 
ного поля: колебания величинь его производньх зависят от флуктуаций малоинтенсивньх 
геофизических факторов (аномальньжх атмосферньмх масс (рис. 2), снеговьїтх масс, уровня 
грунтовьгх вод, лесистости и изменений рельефа (рис. 3) вследствие техногенной деятельнос- 
ти |Двуліт, 1999). Учитьвать зти зффектьт можно путем внесения соответствующих попра- 
вок при решений прямьх задач гравиметрии на исследуемьїх территориях (предполагаєтся, 
что вследствие длительного мониторинга участка приблизительно известно его строение). 

К примеру, влияние притяжения атмосферньїжх масс в центральной сферической зоне 
(сферическая трапеция 192х1?2) вокруг пункта наблюдений на величину гравитационного поля 
определяет формула |Двуліт, 1999) де, - 0.269ДР , где ДР - аномалия атмосферного давле- 
ния; влияние снеговьїх масс в той же центральной зоне равно де, - -21/0, где с - плот- 
ность снега; влияние уровня грунтовьїх вод на вариации Дє учитьваєт формула де 2 2л/0й, 
где п - амплитуда колебания уровня вод; и, наконец, техногенное изменение рельефа - соот- 
ношение де - у - АЛ/Г,, где ДИ - превьшшение между пунктами наблюдений за І год, І, - рас- 
стояние между пунктами. 

Необоснованное упрощение аналитических моделей средьг с целью снижения неоднознач- 
ности истолкования во многих случаях являєтся причиной неверньгх результатов вьтчисле- 
ний геометриий источников, вертикального и латерального распределения плотностньх неод- 
нородностей, особенно в случаях, когда среда вокруг аномального источника далека от пред- 
положений об однородности. Надежное количественноє истолкование динамики масе, на- 
пример, в мониторинге месторождений углеводородов получают при известной геометриий 
тяготеющих тел (по данньм сейсморазведки) и комплексной интерпретации поля силь! тя- 
жести и деформаций рельефа. 

Требование максимального учета имеющейся априорной информации о среде мьг предла- 
гаєм осуществлять двумя путями: 

1. путем построения надлежащих модельньх представлений (звездньєе области известной 
плотности внутри компактньх множеств в банаховом пространстве данньгх); 

2. добавлением в функциональт типа невязки в регуляризирующих алгоритмах некоторьх 
стабилизаторов дифференциального вида, собственньєе функции которьгх совпадают о собст- 
венньми функциями исходньтх операторов. 

Обоснование п. | осуществлено в |Дубовенко, 2002), а п. 2 - в |Черная, 1999). Решениє 
конкретньх обратньх задач методом регуляризации осуществляется по алгоритмам, изложе- 
нньх в |Регуляризирующие..., 1933| и подобньжх им - на основе разработок |Дубовенко, 
2002; Старостенко и др., 2004). 

7. О применениий мониторинга. Повторяемье измерения гравитационного поля применя- 
ют в геодезии, сейсмологиий и вулканологии (комплексное определение предвестников земле- 
трясений и извержений), для мониторинга флюидной динамики водньх бассейнов (рис. 4), 
месторождений и подземньїжх хранилищ углеводородов в процессе их зксплуатации. При 
коррекции методики результать мониторинга можно распространить и на сферу применения 
КПВ, указанную в |Собакадь, 19721. 

При сравнений данньх гравимониторинга месторождения и его аналитической модели 
обнаруживаємое отклонениє наблюденного поля от модели добьтчи является сигналом к 
детальному изучению месторождения методами ГИС и сейсмометрий. Из решения обратной 
задачи гравиметриий по данньм мониторинга можно восстановить общую картину вариаций 
плотности (снижениє Де обозначаєт уменьшение обьема углеводородов вследствиє их до- 
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бьчи и опусканиєс газо-нефтяного контакта, увеличенис Ле - поднятиє уровня пластовьаіх 


вод), но не абсолютньєе значения плотностей. 

Гравитационньй мониторинг миграции СО, в месторождениях углеводородов основан на 
снижений их обьемной плотности при нагнетаниий СО; в месторождение. Так, согласно |Са8- 
регіКоуа, Ноуегяіеп, 2008|, 1-км фронт флюида СО; внутри 20-метрового пласта соли на глу- 
бине 1,9 км в условиях 30905 насьшения СО» и 7096 насьщения раствором соли генерирует на 
поверхности гравианомалию 10 мГал. Решение прямой пространственной задачи с помощью 
| Старостенко и др., 2004| позволяет по данньм микрогравиметрического мониторинга в про- 
цессе нагнетания жидкостей в подземнье хранилища вьіявить их распределение и характер 
миграции по распределению плотностей. 

Большие перспективь имеет решение сопряженной обратной задачи - истолкование вре- 
менньх вариаций гравианомалий, возникающих вследствие изменения положения водонеф- 
тяного контакта либо уровня пластовьїх вод в недрах либо скважинах. Его можно использо- 
вать как недорогой способ гравитационного мониторинга подземньгїх зкосистем мегаполисов 
и решения иньх геозкологических задач (отслеживаний последствий наводнений, оползней, 
динамики загрязнений подземньх бассейнов и т.п.). 

Вв;ражаєм искренню благодарность академику В.И. Старостенко за глубокую дискус- 
сию, содействовавшую значительному углублению и пониманию вопросов данной статьи. 
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